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Abstract 

The new scheme of stochastic quantization is proposed. This quantization procedure is equivalent 
to the deformation of an algebra of observables in the manner of deformation quantization with an 
imaginary deformation parameter (the Planck constant). We apply this method to the models of 
nonrelativistic and relativistic particles interacting with an electromagnetic field. In the first case 
we establish the equivalence of such a quantization to the Fokker-Planck equation with a special 
force. The application of the proposed quantization procedure to the model of a relativistic particle 
results in a relativistic generalization of the Fokker-Planck equation in the coordinate space, which 
in the absence of the electromagnetic field reduces to the relativistic diffusion (heat) equation. The 
stationary probability distribution functions for a stochastically quantized particle diffusing under 
a barrier and a particle in the potential of a harmonic oscillator are derived. 
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I. INTRODUCTION 



There are many different approaches to stochastic quantization and in understanding 
what it is (see for the review jj). In this paper we propose another procedure of stochastic 
quantization, which in some sense generalizes the operator approach to the Fokker-Planck 
equation used in . This new method of quantization gives a stochastic mechanics, 

which is not equivalent to quantum mechanics both in the manner of Nelson's stochastic 
quantization [5J, and the Parisi-Wu stochastic quantization in the fictitious time. Rather we 
interpret the stochastic quantization from the point of view of the deformation quantization 
6|, i.e., as a deformation of an associative algebra of observables (smooth functions over a 
symplectic manifold) with an imaginary deformation parameter as opposed to an ordinary 
quantum mechanics with a real deformation parameter (the Planck constant). This for- 
mulation of stochastic quantization allows us to apply the developed methods of quantum 
mechanics to the stochastic mechanics almost without any changing. 

In this paper we only formulate the general notions of such a stochastic quantization and 
show how it works on simple examples: the models of relativistic and nonrelativistic particles 
interacting with an electromagnetic field. The development of the secondary stochastic 
quantization and its applications to the models with infinite degrees of freedom are left for 
a future work. 

The paper is organized as follows. In the section [III we specify the rules of stochastic 
quantization and introduce basic notions of the proposed stochastic mechanics. In the 
section HTT1 we consider two examples: the stochastically quantized models of a nonrelativistic 
particle in the subsection IIII Al and a relativistic particle in the subsection IIIIBI 

As far as the nonrelativistic case is concerned we find several simple stationary solutions 
to the derived equations of motion: a particle diffuses a potential barrier and a particle in the 
potential of a harmonic oscillator. Here we also obtain the functional integral representation 
for a transition probability and the explicit formula for a first correction to the Newton 
equations due to the diffusion process. Besides we establish that the proposed stochastic 
mechanics can be reproduced by an appropriate Langevin equation. 

In the relativistic case we obtain a Lorentz-invariant generalization of the Fokker-Planck 
equation in the coordinate space, which in the absent of the electromagnetic fields reduces 
to the relativistic diffusion (heat) equation (see for the review [7|). By this example we also 
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show how the basic concepts of the BRST-quantization (see, e.g., 8j) look in the context of 
stochastic mechanics. 

In conclusion we sum up the results of the paper and outline the prospects for possible 
further research. 



II. THE RULES OF STOCHASTIC QUANTIZATION 

In this section we formulate the rules of stochastic quantization and define the main 
concepts of such a stochastic mechanics. 

Let us given a classical system with the Hamilton function H(t, x,p), where x % and pj are 
canonically conjugated with respect to the Poisson bracket positions and momenta 

{x i ,p j } = S i j , i,j = hd, (1) 

where d is a dimension of the configuration space. As in quantum mechanics we associate 
with such a system the Hilbert space of all the square-integrable functions depending on x 
with the standard inner product 

= J d d xr(x)v(x), (2) 

Henceforth unless otherwise stated we consider only real-valued functions in this space. 
In the Hilbert space we define the operators x % and pj such that 

[x\pj] = H8p x l+ = x\ pt = -pj, (3) 

where h is a small positive number and the cross denotes the conjugation with respect to the 
inner product (j2J). Define the Hamiltonian H(t,x,p) by the von Neumann corresponding 
rules 1 

x % — ► x l , pj — > pj. (4) 

The state of the stochastic system is characterized by two vectors and \0) from the 
Hilbert space with the evolution 

h±\1>) = HM, hj t (0\ = -(0\H 1 (5) 



1 We emphasize that contrary to [lj the Hamiltonian H(t,x,p) is not the Fokker-Planck Hamiltonian. 
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and the normalization condition 

(0|V> = 1. (6) 
Define an average of the physical observable T(t,x,p) by the matrix element 

(f) = {0\f(t,x,pM, (7) 

where the operator T(t,x,p) is constructed from T(t,x,p) by the corresponding rules (j4j). 
Then the Heisenberg equations for averages are 

hj t (f) = (d t f+[f,H}}. (8) 

By definition the probability density function is 

p(x) = (0\x)(x\i/>), (9) 

where \x) are the eigenvectors for the position operators corresponding to the eigenvalue x. 
The transition probability from the position x at the time t to x' at the time t' looks like 

G(t',x';t,x) = {0{t')\x'){x'\U tl ,\x)—^- (10) 

(0{t)\x) 

where Ut',t is the evolution operator obeying the equations 

hd t ,u t , tt = HU t ,, u u t ,t = i. (ii) 

The transition probability (flOl) possesses the property of a Markov process 

G(t', x'; t,x) = j d d yG(t', x'; r, y)G(r, y; t, x). (12) 

By the standard means (see, e.g., [9]) we can construct a path integral representation of 
the transition probability ffTOT) . To this end we introduce auxiliary vectors \ip) in the Hilbert 
space such that 

f d d p 

pj\ip) = ipj\ip), (ip'\ip) = 5 d (p-p), / \ip)(ip\ = 1. (13) 



(27rh) c 



In the coordinate representation we have 



% 

(x\ip) = exp{--piX 1 }. (14) 
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Then inserting the unity partition ([13]) into the transition probability (jTUl) we arrive at 

1 



(0(t + dt)\x')(x'\U t+dttt \x) 
(x'\ exp 



dt 
J 



(0(t)\x) 

H(t,x,p + hVlnO(t,x)) + hd t \nO(t,x) J \x) = 
S ^$ eXP {^h [pM^^+^Hit^it + d^M^ + m^O^xit)))] dtj, (15) 
where = x, x{t + <it) = x', x(t) = (x(t + dt) — x(t))/dt, 0(t, x) = (0(t)\x) and 

H(t,x,ip) = (x\H(t,p + JiVlnO(t, x), x)\ip)(ip\x) (16) 

is a gp-symbol of the Hamiltonian H with the momentum p + KVhaO. 

The functional integral representation of the transition probability is obtained by the 
repeatedly use of the property ffT2l) and the formula f|T5|) : 

t) 

x 



<?(*',«';*,») = / n a « n ^ 

J re(t,t') i-6[t,*') 

exp^-- / dr [ Pi (t)x 1 (t) + i(H(t, x(t + dr),ip{r)) + hd T In 0(t,x(t)))] \. (17) 

The property (lT2"j) guarantees that the functional integral representation ([T7fl_does not de- 
pend on what slices the time interval [t,t'] is cut (for more details see, e.g., [10j). 

To conclude this section we formulate the above stochastic mechanics in terms of the 
density operator 

p=mo\. (is) 

From (jnj) and ([6]) it follows that 

hj t p=[H,p), Spp = l. (19) 

The averages are calculated as in quantum mechanics 

(f) = Sp(pT). (20) 

The probability density function p(t,x) is the average of the projector \x)(x\ and obeys the 
evolution law 

hd t p(t,x) = (x\[H,p}\x). (21) 
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As we will see in the next section this equation is nothing but the Fokker-Planck equation. 
Notice that from the definition (j!8p the density operator is idempotent, i.e., 

P 2 = P- (22) 

By analogy with quantum mechanics one can say that such a density operator describes a 
pure state. The transition probability ffTOl) is 

G(t',x';t,x) = Sp(p(t',t)\x')(x'\), p(t,t) = MM (23) 

where p{t',t) obeys the von Neumann equation ffTO]) . 

The formulation of the stochastic mechanics in terms of the density operator reveals 
that from the mathematical point of view the positions x % are not distinguished over the 
momenta Pj as it seems from ([3]). The above stochastic quantization can be considered as 
a formal deformation of the algebra of classical observables in the manner of deformation 
quantization 6] . For a linear symplectic space the Moyal product is 

f(z) * g{z) = J2 ~ ( 2 ^ ■ ■ ■ ^ nbn da 1 ...aJ(z)d bl ... b M, (24) 

n=0 ' ^ ' 

where z = (x,p), a n ,b n = 1,2c?, f(z) and g(z) are the Weil symbols, and u ab is the inverse 
to the symplectic 2-form uj ab . The trace formula for averages is given by 



(T) 



Sp(pT) = J d ^ d ^ d ^/detuj ab p(x,p)T(p,x), (25) 

where p(x,p) and T(p,x) are qp- and pg-symbols of the corresponding operators. For in- 
stance, the gp-symbol of the density operator is 

p(x,ip) = (x\p\ip)(ip\x). (26) 



Thus all the general results regarding deformation quantization of symplectic 11[ and Pois- 



son [12J manifolds, quantization of systems with constraints (see, e.g., |8|) etc. are valid in 
such a stochastic mechanics. 

III. EXAMPLES 

A. Nonrelativistic particle 



In this subsection we consider the stochastic quantization of the model of a nonrelativistic 
particle and in particular establish the one-to-one correspondence of such a quantized model 
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with appropriate Langevin and Fokker-Planck equations. 

According to the general rules expounded in the previous section the Hamiltonian for a 
nonrelativistic particle looks like 2 

(p - A) 2 



H 



2m 



+ A°, 



(27) 



and the evolution equations §5§ in the coordinate representation are 



Hd t ip(t, x) 



(P-A) 

2m 



-A n 



ip(t,x), hd t O(t,x) 



(P + A) 

2m 



-An 



0(t,x), (28) 



where p = — KV and A^(t,x) are gauge fields, which we will call the electromagnetic fields. 
The physical meaning of the fields A^ will be elucidated by the Fokker-Planck equation 
associated with ( J2B1 . 

The equations ( 1281) are invariant under the following gauge transformations 

il>(t,x) -> i>(t,x)e-^> x \ 0(t,x) -> 0(t,x)e^' x \ A^x) -> + <V(t,x). 

(29) 

In particular, these transformations do not change the probability density function. The 
conserved 4-current corresponding to the gauge transformations (|29|) is 



(30) 



f = < 2m" " A) ^ " + A) ° ] 



The system (1281) is Lagrangian with the Hamiltonian action of the form 



S H [0, ip]= dtd d x I hOd t ip - O 



(P ~ A) 2 

2m 



An 



1> 



(31) 



that is the fields if)(t,x) and 0(t, x) are canonically conjugated. 
With the identification 

0{t, x) = e ^ s{t ' x) , x) = p(t, x)e~T, s ^ x \ 



(32) 



the system of evolutionary equations fl28l) becomes 3 
d t p = — div 



h „ V5- A " 
-— Vp+ p 

2m m 



(\7 Q — \ ) 2 h 

d t S-A +^-- L = _ div(VS-A). (33) 

2m 2m 



2 We use the Minkowski metric r/^ = diag(—l, 1, 1, 1) and the system of units in which the velocity of light 
c = 1. The bold face is used for the spacial components of 4- vectors. 



For possible nonlinear generalizations see, e.g., 13 1 
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The first equation in this system is the Fokker-Planck equation, while the second equation 



can be referred to as the quantum Hamilton- Jacob i equation [14]. 

Now it is evident that if one neglects quantum corrections then the initially ^-shaped 
probability density function p(t, x) keeps its own form and propagates as a classical charged 
particle in the electromagnetic fields 4 with particle's momentum VS(t,x) — A(t,x). 

Let us find the evolution of the average position of the stochastically quantized particle. 
The Heisenberg equations (jSJ) for this model are 

m|(x) = (p-A) = (VS-A), m^(x) = (E) + i<(V£- A) x H) + A( rot H ). (34) 

In the case that p(t, x) is sufficiently localized comparing to the characteristic scale of vari- 
ations of the electromagnetic fields the angle brackets can be carried through the electro- 
magnetic fields to obtain a closed system of evolutionary equations on the average position. 
They are simply the Newton equations with the "quantum" correction. 
Notice that the analog of the quantum mechanical uncertainty relation is 

((^) 2 )((PL) 2 )>^ (35) 
where p os = — KV In p 1 / 2 is the osmotic momentum. It is easily proven from the inequality 

J d d x [(fx* - hdi)p 1/2 ] 2 > 0, Vf G R. (36) 



The equipartition law 



151 ] can be discovered from 



l im T{ !^!W dt} = [± h}} = -d, (37) 



where x(t) are the position operators in the Heisenberg representation and T means the 
chronological ordering. 

To reproduce the Fokker-Planck equation associated with the Langevin equation of the 

n 

form (see, e.g., [4J) 

^-x\t) = f%x(t)) + z/(t), (z/(t)> = 0, (^(*y'(0> = hS ij S(t - 0, (38) 
at 

where v % (t) is a Gaussian white noise, one has to solve the system of equations (m = 1) 

VS(t,x)-A(t,x) = f(t,x), A -d t S = ^(f 2 + ftdivf) , (39) 



Such an interpretation for the Langevin equation with a non-conservative force was proposed in 16 1 



s 



with respect to A^(t,x) and S(t,x). Obviously, this system admits a solution. The arbi- 
trariness in the definition of A^t, x) and S(t, x) from the equations (1391) is equivalent to the 
arbitrariness of a gauge. The converse is also true, i.e., for any given solution S(t,x) and 
A^t, x) of the quantum Hamilton- Jacobi equation (1551) we can construct the force f(t,x) 
in the Langevin equation by the formula (1391) . which gives rise to the same probability 
distribution function. The equations (134")) for the average position of the particle in the 
representation fl39|) become 

d d 2 ft 

|(x) = (f), ^(x) = ((d t + (fV))f) + -<Af). (40) 

To gain a better physical insight into the stochastically quantized model of a nonrela- 
tivistic particle we construct the functional integral representation (TlTj) of the transition 
probability. The gp-symbol of the operator appearing in the formula ( fl6|) is 

H{t, x,ip) = — [-p 2 + 2zp( VS - A) - h div( VS - A)] + A . (41) 

Substituting this expression into ffTTj) and integrating over momenta we arrive at 

G(t',x';t,x) 



exp 




Tg(i,t') 

771 n 

-x 2 + (A - VS)± - (A + d T S) - — div(A - VS) 
2 2m 



(42) 



where the functions ^(t,:*;) and S(t,x) obey the quantum Hamilton- Jacobi equation (I3"3"|) 
and are taken at the point (t, x) = (r, x{r+dr)). Now it is obvious that the main contribution 
to the transition probability is made by the paths approximating a classical trajectory. In 
the representation (139]) the transition probability (|42l) reduces to the well known result 

1 T-r d d x(r) 

X 



• (43) 



exp 




(x(r)-f(r,,(r + rfr))) 2 + ^.^^ + ^ 



Usually the force f(t,x) is specified so that the corresponding Fokker-Planck equation 
admits a Boltzmann's type stationary solution. As one can see from the equations fl33l) that 
is the case if VS and A^ are of the order of h or higher, i.e., the momentum and energy of 



the particle are small. For example, the Boltzmann distribution 

p(x) = e~ u ^/Z, 



(44) 



where U(x) is some time-independent potential function measured in terms of the temper- 
ature, is reproduced by the following solution to ( 1331) 



S 



hu 

~2~' 



Ac 



4 



{VU) 2 + AU 



A = 0. 



(45) 



Possibly such "quantum" corrections to the electromagnetic potential naturally arise from 
the stochastic quantization of the electromagnetic fields (we leave a verification of this sup- 
position for future investigations). Nevertheless in a high-energy limit, while the diffusion 
results in small corrections to the dynamics, the gauge fields in the equations ( 1331) can be 
interpreted as the electromagnetic fields. Notice that under this interpretation the equations 
( 133]) are Galilean invariant as opposed to the case, when VS — A is a force. 

To conclude this section we give several simple one-dimensional stationary solutions to 
the equations (|33l) . 

The stationary solutions for A^ = 0. The system of equations ( 1331) is 



hp' = 2S'p, hS" + S 1 



12 



2m£7, 



(46) 



where E is a constant. The solutions are 

,2 



£ = ^>0, 
2m 

£7 = 0, 



E < 0, 



p = Cl e- 2px / h + c 2 e 2px ' h + 2 V ^, 

p=(x — c) 2 /Z, or p = c, 
1 



(47) 



P 



cos 



^2m\E\ 



h 



(x — c) 



In the last case we can take only one hump of the squared cosine function and then continue 
the solution by zero on the residual part of the line. 

To obtain solutions with a finite norm describing a diffusion of particles under a potential 
barrier we just have to join the solutions in f j4Tj) . For a potential barrier of the form 5 



A°(x) = V, x > 0, 



(48) 



5 For brevity, we hereinafter designate only nonvanishing parts of a piecewise function. All the below 
solutions have a continuous first derivative on a whole real line. 
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where V is a positive constant, we have 



1 e^^n+pHlh- 2 )- 1 , x<0; 
P = y { _ _ (49) 

^ I cos 2 {[x — Ip aictan(pl p h 1 )]/^ p }, < x < / p [arctan(p/ p ^, x ) + 7r/2]; 

where < p < (2mV^) 1 / 2 and the characteristic penetration depth 

h 



I, 



(50) 



a/ 2m — p 2 ' 

is of the order of the penetration depth of a quantum mechanical particle (of course, if one 
considers h as the Planck constant). For the potential barrier (|48p there are normalizable 
stationary solutions distinct from (H91) of the form 



I I v x + / O cot(a/Z )) /Iq, x G [— Z cot(a/Zo),0); 

cos 2 [(x — a) /In] r t i \ 

-, x G [0,a + 7T/ /2); 



a G (0, ttZq/2) 



sin 2 (a/Zo) 
For a small potential barrier 



(51) 



A°(x) = V, -1/2 <x< 1/2, I < irl 



o- 



(52) 



we obtain the following stationary solutions 



e 2 P (x+i/2)/ft ; x < _ Z / 2 ; 

cos 2 (x/L) r , . , . , 1 

e -2p(*-»/2)/R j x >//2; 



(x + c) 2 // 2 ,, x G [-c, -Z/2); 
cos 2 (x// ) 



, x G [-1/2,1/2); (53) 



sin 2 [//(2/ )] 
(x-c) 2 // 2 , xG[Z/2,c); 



where p should be determined from the equation p = htan[l/(2l p )]/l p having the unique 
solution and c = l cot[Z/(2/ )] + 1/2. Thus for the barrier of this type the probability to find 
a particle near the barrier is higher than remotely from it. 

The stationary solutions for A = muV/2, A = 0. The system of equations (I33p can be 
rewritten as 



p=Y> f" + (y 2 -e)f = o, 



O(x) 



'(( 



muj \ V 2 



x 



- (54) 



h / J hijj 

Whence from the requirement p'(0) = we have the two types of stationary solutions 

"muj \ 1 / 2 



p = $ ( , -; m/ ) 

^ Z v 4 2' y y 



6 ^ f> T n / 3 'ZfS 3 , 



y 



C- 



v h 



x, (55) 
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FIG. 1: The probability density function for a stochastically quantized particle in the potential of a harmonic 
oscillator. The halves of the first humps normalized on 1/2 are only depicted. The solutions corresponding 
to the first formula in are plotted on the left, while the solutions corresponding to the second formula 
in are plotted on the right. 



where $(a,c;x) is the confluent hypergeometric function (see, e.g., |17j|). As above we can 
take only the part of the solution (1531) defined on the segment between two nearest to the 
minimum of the potential zeros of p(x) and continue it on the residual part of the line by 
zero. It is permissible because p(x) has degenerate zeroes. Then for an arbitrary value of the 
parameter e these distributions are bounded and have finite norms (see Fig. [1]). Otherwise 
the integral of p(x) diverges logarithmically 

It is not difficult to obtain the asymptotic at x — > +oo of a one-dimensional stationary 
solution to (JMD for A = V(x), A = 0: 



P 



cos 2 [ft" 1 J(2mV) 1 / 2 dx] 
ZV 1 / 2 ' 



(56) 



where V — > +oo is assumed. The probability density function p(x) has a finite norm if V(x) 
increases more rapidly than x 2 at both infinities. 



B. Relativistic particle 

In this subsection we stochastically quantize the model of a relativistic particle interacting 
with the electromagnetic fields. As the result we obtain a relativistic generalization of the 
Fokker-Planck equation in the coordinate space. This model also serves as a simple example 
of a model with constraints. 
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The Hamiltonian action for an interacting relativistic particle has the form 6 

S H [x,p, A] = J dr[p^ - \({p - A) 2 + m 2 )], (57) 

where is the electromagnetic potential. The dynamics of the model ( 1ST)) is governed by 
the one constraint of the first kind. 

According to the standard BFV-quantization scheme of the models with constraints of 
the first kind (see, e.g., [8j) we introduce a canonically conjugated ghost pair (c, P) and 
construct the BRST-charge 

Q = c{{p-A) 2 + m 2 ), {Q,Q} = 0. (58) 

The quantum BRST-charge is obtained from the classical one by means of the von Neumann 
corresponding rules (j3j). A graded version of the commutation relations for positions and 
momenta is held. Therefore the quantum BRST-charge is nilpotent but not Hermitian. 
Then the physical state is defined as 

fl|V>) = 0, (O\Cl = 0, [Cl,p] = 0. (59) 

This definition of a physical state respects also the BRST-cohomologies structure, i.e., the 
average over a physical state of a BRST-exact operator vanishes. Explicitly, in the coordinate 
representation we have 

[(hd^ + A M )(n0" + A") + m 2 }ij(x) = 0, [(hd, t - A^Hd" - A") + m 2 ]0(x) = 0. (60) 

When the electromagnetic fields vanish these equations are the Klein-Gordon equations for 
tachyons 7 . 

The action functional for the system of equations (1601) is 

S[0,il>] = J d d xO(x)[(hd fM + A fM )(hd" + A^)+m 2 }i;(x). (61) 

As in the nonrelativistic case the action possesses a gauge invariance under the transforma- 
tions (1291) . The conserved 4-current looks like 

f = i ~ A ^ ~ ^ + A " )0] ' (62) 

6 In this subsection d is a dimension of the space-time and x denotes a set of coordinates on it. 

7 For the interrelation between relativistic random walking models and relativistic wave equations see, for 
instance, 0, 
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where p^ = —Hd^. Making the substitution ([3*2]) into the system (1601) we obtain a Lorentz- 
invariant generalization of the equations fl33l) 



0" 



-~d„p+{d»S-A»)p 



0, {dS - A) 2 + m 2 = -hd^d^S - A^). (63) 



Again the first equation can be called as the relativistic Fokker-Planck equation in the 
coordinate space 8 , while the second equation is the quantum Hamilton- Jacob i equation. In 
the presence of the electromagnetic fields the integral 



/ 



dxp(t,x), (64) 



is not an integral of motion. Analogously to quantum mechanics we can explain it by the 
pair creation. 

In the absence of the electromagnetic fields there is a solution to the quantum Hamilton- 
Jacobi equation (163]) in the form of a "plane wave" 

S = p^x 11 , p 2 = —m 2 . (65) 

Then the relativistic Fokker-Planck equation is rewritten as 

pPdpp = ^Dp. (66) 

That is the celebrated relativistic diffusion (heat) equation (see for the review [3]). It is the 
hyperbolic type differential equation and, consequently, the propagation velocity of small 
fluctuations does not exceed the velocity of light contrary to the nonrelativistic diffusion 
equation. The integral (l64"j) is conserved under an appropriate initial condition. 

Notice that in the same fashion we can quantize the model of a nonrelativistic particle 
in the parameterized form 

S H [x,p, A] = J dr [p^ - \(p Q + H(x,p))} , (67) 

reproducing the results of the previous subsection. 



For the relativistic Fokker-Planck equation in the momentum space see, e.g., 20(. For other approaches 
to a relativistic diffusion see, for example, 21, 22, 2^|. 
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IV. CONCLUDING REMARKS 



There are at least two possible points of view on the results of this paper. 

On the one hand we can consider the proposed quantization scheme from the position of 
deformation quantization. Then we investigate in this paper what happens when the algebra 
of observables is deformed by an imaginary parameter contrary to quantum mechanics with 
the real Planck constant. It would be intriguing if such a deformation results in a stochastic 
mechanics related in some way to real physics. The grounds for these hopes are provided 
by the observation that the obtained stochastic mechanics is closely related to the Langevin 
and Fokker-Planck equations and in the classical limit turns into classical mechanics. 

On the other hand we can regard the proposed quantization procedure as another refor- 
mulation of the Langevin equation. This reformulation treats not only nonrelativistic and 
relativistic models in a uniform manner, but allows us to extend the developed methods of 
quantum mechanics to non-equilibrium statistic physics. 

In both cases the work deserves further research. On this way we can distinguish the 
secondary stochastic quantization and its applications to the models with infinite degrees 
of freedom both in the relativistic and nonrelativistic cases. The most prominent models 
are of course the models of scalar and electromagnetic fields. Then we can attempt to 
attack the model of an incompressible fluid and compare the obtained stochastic model with 
the known one for the fully developed turbulence derived from the Langevin-Navier-Stokes 
equation (see for the review [241]). 
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